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A dispersion relation formalism for the virtual Compton 
scattering (VCS) reaction on the proton is presented, which 
for the first time allows a dispersive evaluation of 4 generalized 
polarizabilities at a four- momentum transfer < 0.5 GeV^. 
The dispersive integrals are calculated using a state-of-the-art 
pion photo- and electroproduction analysis. The dispersion 
formalism provides a new tool to analyze VCS experiments 
above pion threshold, thus increasing the sensitivity to the 
generalized polarizabilities of the nucleon. 
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Over the past years, the virtual Compton scatter- 
ing (VCS) process on the proton, accessed through the 
ep — > epj reaction, has become a powerful and precise 
tool to provide new information on the internal struc- 
ture of the nucleon VCS has been shown to be of 
particular interest not only at low outgoing photon ener- 
gies where one probes nucleon-core excitations and pion- 
cloud contributions to so-called generalized polarizabili- 
ties, but also at high energy and momentum transfers, 
where one is sensitive to a new type of parton distribu- 
tions, generalizing the information obtained from inclu- 
sive deep-inelastic scattering. 

In the low energy regime below pion threshold, the 
outgoing photon in the VCS process plays the role of a 
quasi-constant applied electromagnetic dipole field and, 
through electron scattering, one measures the spatial dis- 
tribution of the nucleon response to this applied field [Q . 
The response is parametrized in terms of 6 generalized 
polarizabilities (CP's) ||^,^, which are functions of the 
square of the virtual photon four- momentum Q^. The 
CP's provide valuable non-perturbative nucleon struc- 
ture information, and have been calculated in different 
approaches In particular, the CP's teach us about 

the interplay between nucleon-core excitations and pion- 
cloud effects. 

The first dedicated VCS experiment has been per- 
formed at MAMI ||l^ and two combinations of CP's have 
been determined at = 0.33 GeV^ Further VCS ex- 
periments are underway at lower at MIT-Bates | |Tl| ] 
and at higher at JLab JT^ . 

At present, VCS experiments at low outgoing photon 
energies are analyzed in terms of a low-energy expansion 
as proposed in , assuming that the non-Born response 



of the system to the quasi-constant electromagnetic field 
of the low energetic photon is proportional to the CP's. 
As the sensitivity of the VCS cross sections to the CP's 
grows with the photon energy, it is advantageous to go to 
higher photon energies, provided one can keep the theo- 
retical uncertainties under control when crossing the pion 
threshold. The situation can be compared to real Comp- 
ton scattering (RCS], for which one uses a dispersion 
relation formalism ||lj,0 to extract the polarizabilities 
at energies above pion threshold, with generally larger ef- 
fects on the observables. The aim of the present work is 
to provide such a dispersion formalism for VCS, as a tool 
to analyze VCS experiments at higher energies in order 
to extract the CP's from data over a larger energy range. 
It will be shown that the same formalism also provides 
for the first time a dispersive evaluation of 4 CP's. 

To calculate the VCS process, we start from the helic- 
ity amplitudes : 
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,{q,X)e:{q\X')u{p\s')M''^u{p,s), (1) 



with e the electric charge, q {q') the four- vectors of the 
virtual (real) photon in the VCS process, and p (p') the 
four-momenta of the initial (final) nucleons respectively. 
The nucleon helicities are denoted by s,s' — ±1/2, and 
M, u are the nucleon spinors. The initial virtual photon 
has helicity A = 0, ±1 and polarization vector e^, whereas 
the final real photon has helicity A' = ±1 and polariza- 
tion vector e^. The VCS process is characterized by 12 
independent helicity amplitudes Tx's';\s- 

The VCS tensor A^^'' in Eq. (0) is then expanded into 
a basis of 12 independent gauge invariant tensors , 



(2) 



as introduced in (starting from the amplitudes of |^ ) . 
The amplitudes Fi in Eq. contain all nucleon struc- 
ture information and are functions of 3 invariants for the 
VCS process : = -q^, v = [s - u)/{4Mn) which is 
odd under s <-> it crossing, and t. The Mandelstam in- 
variants s, t and u for VCS are defined by s = {q+p)'^, 
t = {q — q')^, and u — {q — p')'^, with the constraint 
s + t + u — 2M^ — Q^, and is the nucleon mass. 

Nucleon crossing combined with charge conjugation 
provides the following constraints on the amplitudes Fi 
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1^ at arbitrary virtuality : 

F,{Q\-v,t)^F,{Q^,v,t) (z-l,...,12). (3) 

In a next step, the VCS tensor A^^" at low outgoing 
photon energies is separated into Born [B) and non-Born 
{NB) parts, as described in In the Born process, the 
virtual photon is absorbed on a nucleon and the inter- 
mediate state remains a nucleon, whereas the non-Born 
process contains all nucleon excitations and meson-loop 
contributions, and is parametrized through 6 GP's. With 
the choice of the tensor basis of the resulting non- 
Born amplitudes Fj^^ {i = 1,...,12) are free of all kinc- 
matical singularities and constraints. 

Assuming further analyticity and an appropriate high- 
energy behavior, the non- Born amplitudes Fl^^{Q'^,v, t) 
fulfill unsubtracted dispersion relations (DR's) with re- 
spect to the variable v at fixed t and fixed virtuality : 



Re^;^^(g^^.,^) = -v 

TT 
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with IvcisFi the discontinuities across the s-channel cuts 
of the VCS process. Since pion production is the first in- 
elastic channel, vthr = Jtitt -I- {m^ + t/2 + /2)/{2Mn), 
where m,r denotes the pion mass. 

The unsubtracted DR's of Eq. (||) require that at suffi- 
ciently high energies ^ cxd at fixed t and fixed Q^) the 
amplitudes lmsFi{Q^, v, t) {i — 1,...,12) drop fast enough 
such that the integrals are convergent and the contri- 
butions from the semi-circle at infinity can be neglected. 
The high-energy behavior of the amplitudes Fi is deduced 
from the high-energy behavior of the VCS helicity am- 
plitudes of Eq. (|l]). After some algebra, we obtain the 
following Regge limit for ly oo, at fixed t and : 

Fi, Fg - g2t/"^-(*)-i ^ ^ (5) 



«p(t)-4 ,,QM(t)-3 



(8) 



In Eqs. (Il)-®, we have separately indicated the high 
energy behavior originating from the "pomeron" Q (with 
Regge trajectory ap(t), and ap(0) « 1.08) and from t- 
channel meson-exchange contributions (with Regge tra- 
jectory aAf(O) < 0.5 ). It then follows from Eqs. (^)-(||) 



^We have redefined 4 of the 12 invariant amplitudes of [^] 
by dividing them through v, such that all of them are even 
functions of y. This simplifies the formalism since only one 
type of dispersion integrals needs to be considered then. 

^The pomeron contributes to the sum of helicity conserving 
amplitudes T-^^i.-^^i +T_-^ i ._j i , i.e. gives the same amplitude 
for parallel or antiparallel orientation of the helicities. 



for two amplitudes, Fi and F5, that an unsubtracted dis- 
persion integral as in Eq. (^ does not exist, whereas 
the other 10 amplitudes on the Ihs of Eqs. (^-(||) can 
be evaluated through unsubtracted dispersion integrals. 
This situation is similar as for RCS, where 2 of the 6 in- 
variant amplitudes cannot be evaluated by unsubtracted 
dispersion relations either jlSj . 

To construct the VCS amplitudes Fi and F5 in an un- 
subtracted dispersion framework, one could proceed in 
an analogous way as has been proposed by L'vov [ p^ in 
the case of RCS. The unsubtracted dispersion integrals 
for Fi and F5 are evaluated along the real v-sxis in a fi- 



nite range 



< V < +v„ 



(with i^j, 



1.5 GeV). 



The integral along a semi-circle of finite radius u^^x in 
the complex ^-plane is described by the asymptotic con- 
tribution Ff, which is parametrized by i-channel poles 
(e.g. for — 0, Ff* corresponds to tr-exchange, and F^^ 
corresponds to Tr'^-exchange) . Since the parametrization 
of the asymptotic parts amounts to some phenomenology, 
we will limit ourselves in the present work to those 10 am- 
plitudes Fi (z 7^ 1, 5) for which the integrals in Eq. (^) 
converge, and which do not involve such asymptotic con- 
tributions. A full study of VCS observables within a 
dispersion formalism - requiring, of course, all 12 ampli- 
tudes Fi, and a parametrization of the two asymptotic 
contributions - will be considered in a future work. 

We next consider the non-Born VCS tensor at low en- 
ergy {\q \ 0) but at arbitrary three-momentum \q\ 
of the virtual photon. In this limit, it has been shown 
|^,|| that the non-Born term can be parametrized by 6 
generalized polarizabilitics (GP's), which are functions 
of \q \ and which are denoted by P^p' {\q\). In this 

notation, p {p') refers to the electric (2), magnetic (1) 
or longitudinal (0) nature of the initial (final) photon, L 
(L' = 1) represents the angular momentum of the initial 
(final) photon, and S differentiates between the spin-flip 
{S = 1) and non spin- flip {S = 0) character of the tran- 
sition at the nucleon side. A convenient choice for the 6 
GP's has been proposed in 



p(01,01)0(|^-.|)^ p("-")0(|g|). 



(9) 
(10) 
(11) 



which reduces to the following expressions in the real 
photon limit {\q \ 0) [|| : 

p(oi.oi)0(0) ~ a, P("'")0(0) ~ (12) 
p(oi.oi)i(0) = 0, P("'")i(0) = 0, (13) 
p(OM2)l(o)^^3^ P("'°^)\0) ^ 72 + 74, (14) 

where a (/3) are the electric (magnetic) polarizabilitics, 
and 72, 73, 74 are 3 of the 4 spin polarizabilitics of RCS. 

In terms of the VCS invariants, the limit |(f | ^ 
at flnite \q \ corresponds to v ^ and t — > — at fi- 
nite Q'^ . One can therefore express the GP's in terms of 
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the VCS amplitudes Fi at the point h' = 0, t = 
at finite Q^, for which we introduce the shorthand : 
Fi(Q2) ^ Ff^^ (Q^ ly = 0, t = -Q^) . The relations 
between the CP's and the ^.(Q^) 

can be found in IHI. 
From the high-energy behavior for the VCS invariant 
amplitudes, it follows that one can evaluate the Fi (for 
I 7^ 1, 5) through the unsubtracted DR's 



TT 



I"^^^^-;^— Q^). (15) 



Unsubtracted DR's for the CP's will therefore hold for 
those combinations of CP's that do not depend upon 
the amplitudes Fi and F5. We note however that F5 
can appear in the combination Fc, + AFn, which has the 
high-energy behavior of Eq. leading to a convergent 
integral. Among the 6 CP's, we find the following 4 
combinations that do not depend upon Fi and i<5 : 



p(01,01)0 I 1 p(ll,ll) 
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p(01,01)l 



p(01,12)l 



91 
1 



F2 + {2Fe 



E + Mn 
E 

F12 



1/2 



MnQq 



E + Al 



N 



{iF5 + Fr 
1 



E 

4Fn 



F9) 
1/2 



(16) 



90 

4Ma,Fi2} 



(17) 



p(01,12)l 



V2qa 
{{h + Fr 

V3 



p(ii,ii)i 



E + M 



N 



A En) 



p(ll,02)l 



E 

AMn{2Fq 
E 



1/2 



Mm go 



h)). (18) 



2 6 V / 191 

{go (^"5 + f'y + 4 Fii) + 8 (2 F^ + F^)] , (19) 



where ii^ — p -|- denotes the initial proton cm. 
energy, and go = -^Af ~ E the virtual photon cm. energy 
in the limit |g | = 0. Unfortunately, the 4 combinations 
of CP's of Eqs. (P^)-(P^) can at present not yet be com- 
pared with the data. In particular, the only unpolarized 
experiment [ p^ measured two structure functions which 
cannot be evaluated in an unsubtracted DR formalism, 
as they contain in addition to p(oi>oi)o ;l/2P("'")o of 
Eq. (p^), which is proportional to a -I- /? at = 0, also 
the generalization ol a ^ (3. 

The 4 combinations of CP's on the Ihs of Eqs. (|l6|)- 
( |l9| ) can then be evaluated by unsubtracted DR's, from 
the dispersion integrals of Eq. ( |l5| ) for the FAQ"^). To 
this end, the imaginary parts Im^i^i in Eq. (15) have to 
be calculated by use of unitarity. For the VCS helicity 
amplitudes of Eq. (P (denoted for short by Tf i), the uni- 
tarity equation reads : 



2Im, Tf, = Y,{'^^f5\Px - Pi)Tl^ Tx, 



(20) 



where the sum runs over all possible intermediate states 
X that can be formed. In our present calculation, we 



saturate the dispersion integrals of Eq. ( |15| ) by the domi- 
nant contribution of the ttN intermediate states. For the 
pion photo- and electroproduction helicity amplitudes in 
the range < 0.5 GeV^, we use the phenomenological 
analysis of MAID jl^ , which contains both resonant and 
non-resonant pion production mechanisms. 
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FIG. 1. Imaginary parts of the VCS helicity amplitudes 
of the proton at —t = Q^ = 0.3 GeV'^, evaluated with n'^n 
and 7r°p intermediate states. The labels (A's', As) indicate the 
corresponding helicities. 

In Fig. |l|, we show the imaginary parts of the 12 VCS 
reduced helicity amplitudes n as calculated with ttN 
intermediate states. These reduced amplitudes are ob- 
tained from the full VCS helicity amplitudes of Eq. (|l|) 
by dividing out a common angular factor : 

Tys'Xs = (cos0/2)l'^+^'l (sin 0/2)1^-'^' I n ■ (21) 

In Eq. (pT|), the total helicities in the initial and final 
states are denoted by A = A — s and A' = A' — s' respec- 
tively, and the correspondence between the 12 Ti and the 
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helicity labels is given in Fig. I. The imaginary parts 
were calculated in two different ways, first through the 
helicity amplitudes as expressed in Eq. (]20|). The sum 
over the final states denoted by X contains a phase-space 
integral for the ttN intermediate state, which is then per- 
formed numerically. In a second calculation, the helicity 
amplitudes are first decomposed into a multipole series, 
and the unitarity equation is then implemented for the 
ttN multipoles. It was found that the partial wave expan- 
sion up to orbital angular momentum ^ < 3, is already 
in very good agreement with the numerical integration, 
thus providing a valuable cross-check on the numerical 
calculation. 
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FIG. 2. Dispersion results for 4 of the generalized polariz- 
abilities of the proton (full curves), compared with results of 
0(/) HBChPT [|| (dashed curves) and the 1 near a-model 
(dashed-dotted curves). 

In Fig. H, we show the results for the 4 combina- 
tions of CP's of Eqs. (|l^-(|l|) in the DR formalism, 
and compare them to the results of the 0{p'^) heavy- 
baryon chiral perturbation theory (HBChPT) and 
the linear cr-model M. The ttN contribution to the 
sum p(oi.oi)o + l/2P(ii'")o gives only about 80% of the 
Baldin sum rule because of a non- negligible high- 
energy contribution (of heavier intermediate states) to 
the photoabsorption cross section entering the sum rule, 
which is not estimated here. For the 3 combinations of 
Eqs. (P^)-(P^) of spin polarizabilities, we expect the dis- 
persive estimates with ttN states to provide a rather re- 
liable guidance. By comparing our results with those 
of HBChPT at 0{p^), we note a rather good agree- 
ment for p(oi.i2)i _^ y/^/2p(ii.02)i^ whereas for p(oi.oi)i 
and p(oi.i2)i _ i/(^go)^'^"'"^\ the dispersive results 
drop much faster with Q^. This trend is also seen in 



the relativistic linear cr-model, which takes account of 
some higher orders in the chiral expansion. Since the 
CP's p(iidi)i and p(*'i'i2)i j-eceive non-negligible contri- 
butions from the A and Di^ resonances a complete 
agreement with the ChPT results, dealing a priori only 
with non-resonant excitations, cannot be expected. On 
the other hand, the resonance contribution to p(oi'Oi)i 
is extremely small, and therefore the large deviation of 
the leading order ChPT result indicates the necessity of 
a higher order calculation. 

In conclusion, we have presented a DR formalism for 
VCS and given, for the first time, a dispersive result for 
4 of the CP's of the proton. These evaluations could be 
used to check the convergence of the chiral expansions in 
ChPT calculations. It also provides a new tool to analyze 
VCS experiments at higher energies (above pion thresh- 
old) where there is an increased sensitivity to the CP's 
of the proton. In the future we will study further details 
of such an analysis, including possible parametrizations 
of the 2 non-convergent dispersion integrals. 
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